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ABSTRACT 

Using GALFORM, a semi-analytical model of galaxy formation in the A cold dark matter 
cosmology, we study the rest-frame ultraviolet (UV) colours of Lyman-break galaxies 
(LBGs) in the redshift range 2.5 < z < 10. As the impact of dust on UV luminosity 
can be dramatic, our model includes a self-consistent computation of dust attenuation 
based on a radiative transfer model. We find that intrinsically brighter galaxies suffer 
stronger dust attenuation than fainter ones, though the relation has a large scatter. The 
model predicts galaxies with UV colours consistent with the colour selection regions 
designed to select LBGs in observational surveys. We find that the drop-out technique 
that selects LBGs based on two rest-frame UV colours is robust and effective, selecting 
more than 70 per cent of UV bright galaxies at a given redshift. We investigate the 
impact on the predicted UV colours of varying selected model parameters. We find 
that the UV colours are most sensitive to the modelling of dust attenuation and, 
in particular, to the extinction curve used in the radiative transfer calculation. If 
we assume a Milky Way dust extinction curve, the predicted UV continuum slopes 
are, in general, bluer than observed. However, we find that the opposite is true when 
using the Small Magellanic Cloud dust extinction curve. This demonstrates the strong 
dependence of UV colours on dust properties and highlights the inadequacy of using 
the UV continuum slope as a tracer of dust attenuation without any further knowledge 
of the galaxy inclination or dust characteristics in high redshift galaxies. 

Key words: galaxies: evolution, galaxies: formation, galaxies: high- redshift. 



1 INTRODUCTION 

The recent availability of near infrared (IR) imaging on 
the Wide Field Camera 3 (WFC3) on board the Hubble 
Space Telescope (HST) , extending the wavelength coverage 
to span 2000A to 17000A, has dramatically increased the 
number of candidate high redshift galaxies with z > 2.5 (e . < 
Hathi et al.ll20ld: lOesch et alJl201Cl: IWilkins et alj|2011al li 



McLure et al.ll201ll : iBouwens et alj|2012h . The galaxy red- 
shift "record" has been tentatively pushed up to z = 10, with 
3 candidates awaiting spec troscopic confirmation l|Yan et alj 
l20ld : lBouwens et alJl2oTlh . 

Lyman-break galaxies are selected by targeting spec- 
tral features of star-forming galaxies at z > 2.5, in particu- 
lar the Lyman break, a combination of both the Lyman- 
limit break (at 912A) due to neutral hydrogen in each 
galaxy and the Lyman-a-break (at 1216A) due to absorption 
by in tervening neu t ral ga s clouds, and the UV continuum 
slope. ISteidel et all (11996) proved the usefulness of this tech- 
nique by spectroscopically confirming drop-out candidates 
isolated photometrically at z ~ 3. Since then, this approach 



has been adapted to find higher reds hift galaxies using 
both ground and space based telescopes (ISteidel et alJll999l; 
Bouwens et all 120031 : IShimasaku et al.l 120051 : lYoshida et al.l 



20061) . 



At these high redshifts, we are witnessing the rapid evo- 
lution of very young galaxies, when the Universe was less 
than a quarter of its current age. Studies of galaxies se- 
lected using the drop-out technique have been fundamental 
in unveiling the global star formation h istory of the Universe 
|Madau et al.lll996l : ISteidel et al.lll999l ). Moreover, observa- 
tions at redshifts above z — 6 could help us to understand 
the end of the reionisa tion of the intergalactic medium (e.g. 
iRobertson et al.ll201ol ). However, the interpretation of the 
observational data is likely to be significa ntly affected by 
the a ttenuation of starlight due to dust (e.g. lChapman et al.l 
l2005h . 

Currently, there is some con troversy over the lev el of 
dust attenuation inferred in LBGs. iMcLure et aLl(|201ll ) and 
iFinkelstein et all (|201lf ) estimated the far UV attenuation 
using spectral energy distribution (SED) fitting for LBGs 
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brighter than M A s(1500A) - 51ogh = -17. Both studies 
used template spectr a including dust attenuation follow- 
ing the Calzetti law. iFinkelstein et all inferred an attenu- 
ation at 1500 A, Aisoo, of ~ 1.3 magnitudes at z ~ 4 
and Ai5oo < 0.25 at z ~ 7. In contrast, for galaxies with 
(z) — 6.5, iMcLure et al.l fo und (A1500) ~ 0.4, a value above 
the upper limit found by IFinkelstein et all In many ob- 
servational studies the dust attenuation is inferred from 
the UV continuum s lope estimated from a single colour. 
iBouwens et al.l (|201 lh measured an average UV continuum 
slope of -3 for galaxies at z = 7. However, this value was 
measured to be ~ —2.2 when more d ata were collected by 
the HST WFC3 |Bouwens et alj|2012h . This illustrates how 
the scarcity of high redshift data can bias the estimation of 
dust attenuation. 

Dust can absorb and scatter photons either into or out 
of the line of sight to the observer. The outcome of these pro- 
cesses depends strongly on the geometry of dust and stars 
in galaxies and a full radiative transfer model is required 
for a realistic calcula t ion of the effect of dust on s tar light 
|Fontanot et alj|2009l ; iFontanot fc SomerviHell201ll ). Obser- 
vationally, the attenuation by dust of the UV light emit- 
ted by LBGs is usually estimated f r om the observ e d UV 
continuum slope (|Meurer et al.ll 19991 ). ICalzetti et ail (|200Ch 
proposed a universal shape for the dependence of the at- 
tenuation on wavelength, such that the UV light is the most 
affected by dust. This so called "Calzetti law" is based on ob- 
servations of 39 nearby UV-bright starburst galaxies. How- 
ever, the Calzetti law provides a poor fit to the UV colours 
of larger samples of nearby galaxies that have a range of 
star formation r ates (jConrov et alj|2010l : iBuat et al. I l2011bl ; 
IWild et al.ll201lh . Other observational studies also find that, 
at least for part of the UV-selected galaxies at 1 < z < 4, 
the attenuation is not well described by the Calzetti law 
(iNoll et al.ll2009l ; IBuat et alj2011al , l2012j ; lReddv et al.ll2012l ; 
iLee et al.ll2012h . 

In this paper we use a semi-analytical approach to 
model the formation and evolution of galaxies with a realis- 
tic treatment of dust. In particular, we study the rest-frame 
UV colours of LBGs within a cold dark matter universe us- 
ingGALFORM, the se mi-analytical model developed initially 
bv lCole et ail (|2000t ). 

Previous th eoretical studies o f LBGs include both 
semi-analytica l (IBaugh et alj [l99Sl; iGuo fc White! 120091; 
Lo Faro etall 120091 ; ILacev et al.l 1201 ll ; ISomerville et all 



2012T) and gas-d ynamical simulations of galaxy formation 



..Nagamind 120021; iFinlator et al.ll2q06l; iForero-Romero et al.l 
l2010l ; ICenll201ll ; |jaacks et al.ll20T2l : iDaval fc Ferrarall2012f ). 

This latter type of simulation has the drawback that, on 
the whole, they do not predict a present-day galaxy popula- 
tion consistent with observations, unlike the semi-analytical 
models. Furthermore, most of these models treat the at- 
tenuation b y dust by using eith e r the Calzetti law o r the 
slab model. iGuo fc White! (|2009h . lLo Faro et al.l (|2009h and 
ISomerville et alj (|2012F investigated the effects on the in- 
ferred luminosity functions and other properties of apply- 
ing observatio nal LBG colour selections to model galaxies. 
IGuo fc White! u sed a phenome nological model for dust at- 
tenuation, while ILo Faro et al.l used a physical model simi- 
lar to that in the present p aper, but restrict ed their analy- 
sis to LBGs at 4 < 2 < 6. ISomerville et al.l studied galax- 
ies at < z < 5, calculating the dust attenuation using a 



slab model with two dust components, diffuse and molecular 
clouds. They found that in order to reproduce the observa- 
tions, they needed to let the effective optical depth normal- 
isation evolve with redshift. 

In this study, we use the version of GALFORM developed 
bv lBaugh etalli|2005h model, but with the theoretical treat - 
ment of dust attenuation described in I Lacev et all i|201lf ). 
in place of the GRASIL (|Silva et al.lll99a ) model, which was 
used in the original paper. GRASIL is a spectrophotometric 
code which is quite expensive computationally and, thus, 
difficult to run for large samp les of galaxies. The parame- 
ters of the iBaugh et al.l (|2005r ) model were set in order to 
match observations of sub-mm selected galaxies, with a me- 
dian redshift z ~ 2, and the LBG luminosity function at 
z ~ 3 within a single framework, at the same time as repro- 
ducing the pro perties of loc al galaxies. ILacev et al.l (|201lD 
found that the IBaugh et al.l model predicts UV luminosity 
functions that agree with observations up to 2 = 10, the 
highest redshift for which measurements are currently avail- 
able. This agreement is notable as most of these observations 
were not available at the time the model was proposed. A 
typical attenuation in the UV of about two magnitudes was 
predicted for bright galaxies, posing the questions of whether 
or not a model with such strong attenuation requires galax- 
ies with unrealistic dust contents, and if it could produce 
UV colours consistent with those observed. 

Here, we show that the Baugh et al. model produces 
galaxies with the right colours to be selected as LBGs over 
the redshift interval 3 < 2 < 10. We also investigate the at- 
tenuation by dust in the model. Intrinsically bright galaxies 
are predicted to be heavily extincted in the model. Neverthe- 
less, the UV continuum slope inferred from the attenuated 
photometry of model galaxies is not necessarily redder than 
observed. In fact, our calculations show that the UV con- 
tinuum slope is a poor indicator of dust attenuation, as our 
results are extremely sensitive to the choice of extinction 
curve used as the input to the attenuation of starlight cal- 
culations. For our standard choice, a Milky Way (MW) ex- 
tinction curve, galaxies get bluer when they are attenuated, 
which, as we explain in Section JB] is due to the presence of 
a bump at 2 175 A in the extinction curve. The discrepancy 
between the model predictions and the observations for the 
slope of the attenuated UV continuum slope is smaller than 
the difference in the predictions on replacing the MW ex- 
tinction curve by an SMC-like extinction curve. 

In section &\ w e summarise the main features of the 
IBaugh et al.l (|2005l ) model. The drop-out technique is in- 
troduced in section Ej3j Section i|4] presents the predicted 
dust masses and attenuation in the UV for bright galaxies 
at different redshifts. In sections ij5] and ij6] we explore the 
predicted UV colours and UV continuum slope for galaxies 
with 2 > 2.5. The conclusions may be found in SJ7] 

All the magnitudes quoted are on the AB system, unless 
stated otherwise. 



2 GALAXY FORMATION MODEL 

We predict the UV luminosities and colours of galaxies in 
a A Cold Dark Matter (ACDM) universe using the GAL- 
FORM semi- analytical galaxy formation model initially de- 
veloped by ICole et all |2000l ). Semi-analytical models use 



© 0000 RAS, MNRAS 000, 000-000 



LBGs UV colours and attenuation 3 



simple, physically motivated equations to follow the fate of 
baryons in a universe in which structur e grows hiera r chically 
throu gh gravitational instability (see iBaughl 120061 ; iBensonl 
|2010| , for an overview of hierarchical galaxy formation mod- 
els). 

GALFORM models the main processes which shape the 
formation and evolution of galaxies. These include: (i) the 
collapse and merging of dark matter haloes; (ii) the shock- 
heating and radiative cooling of gas inside dark matter 
haloes, leading to the formation of galaxy discs; (iii) qui- 
escent star formation in galaxy discs; (iv) feedback from su- 
pernovae (SNe) and from photoionization of the intergalac- 
tic medium (IGM) ; (v) chemical enrichment of the stars and 
gas; (vi) galaxy mergers driven by dynamical friction within 
common dark matter haloes, leading to the formation of 
stellar spheroids, which also may trigger bursts of star for- 
mation. The end product of the calculations is a prediction 
for the number and properties of galaxies that reside within 
dark matter haloes of different masses. 

We focus our attention on the iBaugh et all (120051 ) 
model. Some of the key features of this model are (i) a time- 
scale for quiescent star formation th at varies simply a s a 
power of the disc circular velocity fsee lLagos et al.ll201ll . for 
a study of different star formation laws in quiescent galax- 
ies), (ii) bursts of star formation are triggered only by galaxy 
mergers, (iii) a lKennicuttl (|l983l ) initial mass function (IMF) 
is adopted in quiescent star formation in galactic discs, while 
in starbursts a top-heavy IMF is assumed, with slope x = 0, 
(iv) the inclusion of SN e feedback with the possible occur- 
rence of superwinds (see [Benson et al.ll2003l . for a discussion 
of the effect that feedback has on the luminosity function of 
galaxies), and (v) feedback from the ionisation of the IGM 
is approximated by a simple model in which gas cooling is 
completely suppressed in haloes with circular velocities less 
than 30km s _1 at redshifts z < 10 (|Lacev et al.l 1201 ll ). In 
sections H2.1l and i\2.2\ we give a detailed account of how the 
attenuation due to dust and the IGM are implemented in 
this model. The parameters of this model were fixed with 
reference to a subset of the av ailable observa tions of galax- 
ies, mostly at low redshift. The lBaugh et all model uses the 
canonical (ACDM) parameters: matter density, fio = 0.3, 
cosmological constant, Ao = 0.7, baryon density, Qt = 0.04, 
a normalisation of density fluctuations given by <rg = 0.93 
and a Hubble constant today of Ho = 100ft kms Mpc , 
with h = 0.7. The model employs merger trees generated us- 
ing a Monte Carlo algorithm. Here w e use the algorithm in - 
troduced bv lParkinson et al.l (|200Sl ). iGroizalez et all (|201ll ) 
reported that using this particular approach does not change 
an y of the orig i nal m o del predictions. We refe r the reader 
tolBaugh et~atl (|2005t ). lLacev et al.l (|200Sl ) and lLacev et all 
|201ll ) for a full description of this model. 

In addi tion to reproducing local galaxy data, the 
IBaugh et al.l model matches the number and redshift dis- 
tribution of galaxies det ected by the i r emi ss ion at sub- 
millim etre wavelengths (|Baugh et al.l 120051 : lLacev et alj 
|200S| ), the rest-frame UV 
break galaxies up to z ~ 



predicts that most present day galaxies with masses close to 
that of the MW had a LBG progenitor at 3 < z < 4. 

No parameters have been tuned for the study presented 
here. We do vary some of the parameters in the model in SJS] 
and ^6] to illustrate the sensitivity of the model predictions 
to our choices. 



2.1 Dust attenuation 



Following the approach p resented in lCole et al.l (|2000l ) and 
lLacev et al.l (|201ll . |2012| ). the attenuation of starlight by 
dust is modelled in a physically self-consistent way, based on 
the results of a radiative transfer calculation for a realistic 
geometry in which stars are distributed in a disk plus bulge. 
This geometry is the same as that assumed in m ore detailed 
comp utations carried out in the GRASIL code (|Silva et alj 
1 1998T ) . An important feature of the method is that the dust 
attenuation varies self-consistently with other galaxy prop- 
erties such as size, gas mass, and metallicity. 

Dust can absorb photons and scatter them either into 
or out of the line of sight to the observer. Here we refer to 
dust "attenuation" as the average effect that dust has on 
starlight once the geometric configuration of dust and stars 
is taken into account. Dust attenuation includes the effects of 
both light absorption and scattering, which strongly depend 
on how dust is distributed with respect to the stars in a 
galaxy. We can express the dust attenuation, a^, as the ratio 
between the attenuated and unattenuated luminosity of a 
galaxy: 



L\ (attenuated by dust) 
L\ (unattenuated by dust) 



(1) 



(|Lacev et al.l 12011 


) and the 


vman-a emitters ( 


Orsi et al.l 



The effective optical depth at a given wavelength can be 
related directly to the dust attenuation at that wavelength: 
Toff, a = ~hi(aA). We can also express the attenuation in 
terms of magnitudes, Aa, and then r c ff.A = AA/(2.51ogej]]. 

Usually the term dust "extinction" refers to the attenu- 
ation of the light from a point source placed behind a screen 
of dust. Thus, the "extinction" is independent of the geom- 
etry of the system. 

In GALFORM, dust is assumed to be in two components: 
dense molecular clouds embedded in a diffuse component. 
Stars inside molecular clouds have their light attenuated 
both by the clouds and by the diffuse medium, while the 
remaining stars are attenuated only by the diffuse medium. 



2.1.1 Diffuse component 

The attenuation of starlight by the diffuse component of 
dust is calcu l ated u sing the radiative transfer models of 
iFerrara et ail |l999l ). These models assume that stars are 
distributed in two components: an exponential disk and a 
bulge; and that dust is smoothly distributed in the same 
plane as the stell ar disk. Given an exti nction law, by default 
the MW one, the lFerrara et all (| 19991 ) models give dust at- 
tenuation factors separately for the disc and spheroid light, 
as a function of wavelength, galaxy inclination, ratio of bulge 



to explore the emiss ivity of ionising photons at z > 6. 
IGonzalez et~aT1 (|2012l ) showed that the IBaugh et alj model 



1 Throughout all the paper we use the notation log to designate 
the logarithm in base 10. 
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Figure 1. The optical depth normalised by its value in the 
V-band as a function of wavelength. We plot the extinction 
curves from the MW (solid line) and SMC (dashed line) used by 
iFerrara et al . (1999) to estimate the attenuation curves of galax- 
ies, and the attenuation curve proposed bv lCalzetti et al] ||2000| ) 
to describe observations of nearby star forming galaxies (dotted 
line). 



to disc radial dust scalelength, r e /hii, ratio of dust to stel- 
lar vertical scaleheights, h z ^ uat /h ZtStara , and Ty . Here Ty is 
the extinction optical depth in the V-band looking through 
the centre of a face-on galaxy. In GALFORM, the diffuse com- 
ponent of dust is assumed to follow the distribution of stars 
in a disk or in a starburst component. By default, the dust 
is assumed to have the same radial and vertical scalelengths 
as the stars. The rest of the above properties are predicted 
by GALFORM. 

In order to calculate the extinction optical depth in the 
V-band, rv , the extinction Ay along the line-of-sight is de- 
termined from the neutral hydrogen column density, Nu, 
and the cold gas metallicity, Z co \a'- 



A v = 



Ay 



(2) 

^0 / V V H/ Q 

The above extinction is normalised in such a way that 
gas with solar metallicity, Zq — 0.02, has the lo- 
cal interste llar medium 
iGuiderdoni 



(ISM) dus t-to-gas 
Rocca- Volmerangeil 198?1 ) 



ratio (e.g. 
The neutral hydro- 



gen column density is related to the projected central surface 
gas mass in the disk, M co id/(27r/ifj), and thus: 

' M coM pc 2 1 



0.043 



2irh% Mq 



0.02 



(3) 



The extinction optical depth, ry- , is obtained in the same 
way independently of which extinction curve is being used 
afterwards for the calculation of the attenuation. 

Fig- [U shows the normalised optical depth in the UV to 



near IR wavelength range, f or the extinction cur ves from 
the MW and SMC used by IFerrara et all (|l999l ) to esti- 
mate the attenuation curves of galaxies. For comparison, 
Fig.[T]also displays the attenuation curve empirically derived 
fro m observations of nea rby UV-bright, starburst galaxies 
by ICalzetti et~aH (|2000l ). For wavelengths A > 0.35/xm, 
the three curves are roughly consistent with each other. 
However, at smaller wavelengths the SMC extinction curve 
displays much higher extinction than both the MW and 
Calzetti curves. The MW extinction curve presents a bump 
at A ~ 2175A that is not seen in the other two curves. The 
presence of this bump is related to the existence of small 
carbonaceous particles t hat could be destroyed b y the am- 
bient UV radiation field |Fischera fc Dopita[|201ll ). Though 



less pronounced, such a bump has been o bserved in the 
Large Magellanic Cloud (|Gordon et alJlioOSf ) and in nearby 
star-forming gal axies (over 23000 galaxies at z ~ 0.07, 
IWild et alll201ll ). However, the bump at A ~ 2175A was 
not det ected in the sample of 39 starburst galaxies anal- 
ysed bvlCalzetti et all |l994 12OO0I ). It is worth noting here 
that lViih et al.l(|200St ). using a radiative transfer calculation, 
found the Calzetti attenuation law curve in the particular 
case of assuming ry = 1.5 and a clumpy dusty shell with a 
SMC extinction curve. 



2.1.2 Molecular clouds 

In the model, stars are assumed to form in molecular clouds 
from which they will eventually leak out. The attenuation 
of starlight by dust in molecular clouds is calculated analyt- 
ically, assuming new born stars to be a point source in the 
centre of a spherical cloud with uniform density. 

The attenuation by clouds depends both on the fraction 
of light produced within the cloud and the extinction optical 
depth. The light produced within molecular clouds depends 
on the star formation and chemical enrichment history, the 
IMF and the wavelength. The attenuation by clouds mostly 
affects the UV light, which is mainly produced by young 
stars with a short lifetime. Thus, the fraction of light pro- 
duced within a cloud is obtained by assuming the recent star 
formation rate to be either constant in the case of quiescent 
disks or to decay exponentially for bursts of star formation. 

The extinction optical depth is calculated from the col- 
umn density of gas through a cloud, assuming the same ex- 
tinction law as for the diffuse component of dust. 

As expected, the attenuation calculated in this way de- 
pends on wavelength, the metallicity of the gas and the sur- 
face density of the clouds, which, in turn, depends on the 
cloud mass and radius. For further details see Lacey et al. 
(2012, in prep). 

The values of the parameters affecting the dust attenu- 
ation cal culation in mo lecular clouds are set to be the same 
as in the iBaugh et al.l model. The time for stars to escape 
from their birth molecular clouds is set to t csc = 1 Myr 
and the fraction of dust in molecular clouds to f c i ou d = 
0.25. The value of f c i ou d is consistent, within a factor of 2, 
with estimations bas ed on observations of nearby galaxies 
|Granato et al.l [20001 ). The value of t eS c was chosen in or- 
der to reproduce the luminosity function of observed LBGs 
at z = 3. Moreover, the chosen value for t esc is consistent 
with the escape time inferred from observations of stars in 
the MW jLeisawitz fc Hauserlll988l ). By default, the mass 
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of clouds is set to be M c ioud = iu iviq 
to r c ioud =16 pc , to match the choices made in GRASIL 
l|Silva et al.l 1 1998ft . that, in turn, were made to match ob- 
servations of giant molecular clouds in the MW and to fit 
the spectral energy distributions (SED) of nearby galax- 
ies. With these parameters the far-UV light of stars within 
molecular clouds is almost completely extincted. Thus, the 
net attenuation by clouds is insensitive to the exact value of 
the surface density of a single cloud, which is proportional 
to Mdoud /rcioud! and it will only depend weakly on the as- 
sumed shape of the extinction curve used in the calculation 
of dust attenuation in molecular clouds. 



2.1.3 Dust to metals ratios 

In GALFORM, a fixed fraction of the total heavy elements 
in the gas component is assumed to be locked up in dust 
grains. The dust to metal ratio is thus also fixed in the 
model. The model uses an instantaneous recycling approx- 
imation. This approximation is reasonable as most metals 
are produced by Type II SNe. Even at high redshifts, an ap- 
preciable difference i n returned mass fractio n only persists 
for less than 10 7,5 yr jNagashima et al]|2005ft . No attempt is 
made to model in detail the formation and destruction of 
dust grains. 

In principle, the dust to metal ratio depends on the 
dust production mechanism, which in turn depends directly 
on the star formation history of a galaxy. A consensus on the 
dust production mechanism at different redshifts has not yet 
been reached and thus, assuming a fixed dust to metal ratio 
is a reasonable approximation. At low redshift, the bulk of 
the dust is thought to be genera ted by stars in the as ymp- 
totic giant branch (AGB) phase. IValiante et al.l (|2009h esti- 
mated that even at z ~ 6.4 between 50 and 80 per cent of 
the total dust mass could be generated by AGB stars. This 
percentage depends strongly on the assumed star formation 
history. AGB stars need to reach on average an age of 400 
Myr before they produce dust. Thus, at higher redshifts, 
2 > 6.5, there is not enough time for typical AGB stars to 
generate the dust that has been observed in some systems 
(Dwek et al1l2007l ; [Michalowski et al.ll2010ft . Although SNe 
can gen erate large amounts of du st they also destroy it in 
shocks |Bianchi fc Schneider! 12007ft . The most plausible ex- 
planation for the observed dust at z > 6.5 appears to be the 
dust gr owth in the ISM generated within dense molecular 
clouds jMichalowski et~al]|2010l ; iDwek fc Cherchnefj [201ll ; 
IValiante et alj|201lft . 



2.2 IGM attenuation 

The neutral hydrogen in the IGM absorbs and scatters pho- 
tons from both the Lyman series (bluewards of 1216A) and 
the Lyman continuum^. Within GALFORM this effect can be 
impl emented follow in g the attenuati on prescriptions of ei- 
ther [Madau| | | 1995ft or IMeiksinl (|2006ft . The predicted atten- 
uation at 1000A due to the IGM increases from about 0.5 



2 The absorption due to the Helium present in the IGM is, in 
general, negligible compared with that by neutral hydrogen at 
the wavelengths of interest here. 



magnitudes at 2 ~ 4 to 2 magnitudes at 2 ~ 6, with either 

pres cription . 

iMadaul (| 1995ft used an empir ical approach to model the 
IGM attenuation, while IMeiksinl based part of the prescrip- 
tion on the results of a ACDM simulation with similar cos- 
mological parameters to those used in the B05 model. Nu- 
merical simulations have proved to be successful at repro- 
ducing different observational properties of absorbing sys- 
tems. 

The effective optical depth through the IGM depends on 
the number of absorbers, basically neutral hydrogen clouds, 
along the line of sight. Each absorber is characterised by its 
redshift, its neutral hydrogen column density, Nm, and its 
Doppler parameter, which is related to the temperature and 
kinematics of the cloud. The absorption of the Lyman con- 
tinuum is mainly caused by Lyman limit systems (LLSs), 
with iVm ~ 10 17 cm -2 and scattering by the Lyman a for- 
est (LAF), optically thin systems at the Lyman edge, with 
JVhi ~ 10 13 c m~ 2 

IMadaul and IMeiksinl assumed different evoluti on with 
redshift for the number of LAF and LLSs. While IMadaul 
consider ed that the LAF followed a Poisson distribution, 
IMeiksinl obtained the LAF distribution from a N-body cos- 
mological simulation. 

An empirica l evo lution of LLSs with redshift is assumed 
by both IMadaul and IMeiksinl . However, only the latter is 
consistent with current observations in the r edshift range 
0-3 < 2 < 4 (see Fig. 3 in llnoue fc Iwatall2008ft . 

Ano ther important difference between the IMadaul and 
IMeiksinl modelling of the intergalactic attenuation is that 
in the fo rmer a fixed Doppler parameter is assumed, while 
IMeiksinl varies this value based on simulations , a change 
supp orted by several observa tional studies (e.g. iKim et all 
Il997ft . Moreover, the IMadaul model almost always predicts 
a smaller transmission than empirical models based on th e 
most recent observations (see Fig. 6 in llnoue fc Iwatall2008ft . 

In light of the above discussion, the results p resented 
in thi s paper are calculated by default using the IMeiksinl 
l|2006ft equations. 



3 LYMAN-BREAK SELECTION FOR 
GALAXIES AT HIGH REDSHIFTS 

Observationally, the most extensively used method for se- 
lecting star-forming galaxies at high redshifts is the Lyman 
drop-out technique (fo r details on other methods see the 
review bv lDunlodr2012ft . This technique uses near-UV, opti- 
cal and near-IR filters (depending on the target redshift) 
to sample the redshifted Lyman-break, thereby selecting 
star-forming galaxies, the so called Lyman Break Galaxies 
(LBGs). This method has been commonly applied to select 
galaxies with 2 > 2.5, since their Lyman-break is shifted be- 
yond the observed frame UV range that is difficult to observe 
with ground-based telescopes. 

The drop-out technique typically uses 2 colours, one 
from filters bracketing the Lyman-break and the second 
quantifying the UV continuum slope. The rest-frame UV 
emission, 300 < A(A)< 3000, is largely due to the presence 
of young and massive OB stars in galaxies. Thus, the rest- 
frame UV emission of galaxies that are not actively form- 
ing stars is, in general, negligible in comparison with that in 
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Figure 2. An illustration of the drop-out technique used to select galaxies at 3 < z < 3.5. Left •panel: The intrinsic sp ectrum of a 
star-forming galaxy at z = 3.2, computed with the PEGASB.2 population synthesis model jFioc fc Rocca-Volm crangc 1993) (black line, 
no attenuation is included) plotted together with th e response functions of the U n , G and TZ filters (grey lines). Central panel: The 
(Un-G) versus (G-1Z) colour-colour plot proposed by ISteidel et al.l il995f) for selecting galaxies at 3 < z < 3.5 (grey region) with the 
predicted density contours for galaxies at different rcdshifts, as indicated in the legend. The density contours enclose 80 per cent of the 
galaxies brighter than TZ < 25.5 and G < 27.29 at each redshift. Right panel: The predicted completeness (top panel) and luminosity 
function (main panel) at z = 3.2 of all galaxies (black line) and galaxies with TZ < 25.5 and G < 27.29 and the colour cuts shown in the 
central panel as a grey region (red lines). The redshift dis tribution of this last subsam ple is shown in the inset. For comparison, we show 
as green symbols the observed luminosity function from ISawicki fc Thompsonl lj200rj l (triangles, 1700A), and iReddv fc Steidell {2009) 
(circles, 1700A). 



star- forming galaxies. The Lyman-break is the result of both 
the Lyman-limit break (at 912A) due to neutral hydrogen in 
the atmospheres of massive stars and the Lyman-Q-break (at 
1216A) due to absorption by intervening neutral gas clouds 
(e.g. the Lyman-a? forest). 

The drop-out techn ique was first proposed by 
ISteidel fc Hamilton! j 19921 ). using the U n , G and TZ niters 
shown in Fig. O The grey area in the central panel of Fig . 
[2] shows the selection criteria used by ISteidel et all (| 19951 ) 
to select galaxies at 3 < z < 3.5. The left panel in Fig. [2] 
shows the synthetic spectrum of a star-forming galaxy at 
z = 3.2, comput ed with the PEGASE.2 stellar p opulation 
synthesis model (|Fioc fc Rocca-Volmerange] Il997l ) without 
including the attenuation by either dust or the IGM. This 
spectrum shows a very pronounced break at 912A in the rest 
frame, the Lyman-limit break. Bluewards of this rest-frame 
wavelength, the galaxy is practically not emitting any light. 
Thus, by measuring the (U„-G) colour of a galaxy it is possi- 
ble to select candidates at 3 < z < 3.5, since in this redshift 
range the U n band is bluewards of the Lyman-break and the 
G band redwards. The (G-TZ) colour samples the rest-frame 
UV continuum slope and is used in order to minimise the 
number of lower redshift interlopers. 

The central panel of Fig. [2] shows, using density con- 
tours, the location in the colour-colour plane of the model 
galaxies at different redshifts. It is clear from this plot that 
a selection based only on the colour sampling the Lyman- 
break (i.e. (U n -G) in this example) will contain a substan- 
tial fraction of galaxies at z ~ 0. Old and passively evolving 
galaxies at z ~ present a deep enough break at 4000A to 
be mistaken for the Lyman-break of star-forming galaxies at 
higher redshifts. However, the slopes of the spectral contin- 
uum redwards of these breaks are very different. The rest- 
frame optical continuum of old galaxies at z ~ is quite red, 



while the rest-frame intrinsic UV continuum of star-forming 
galaxies is bluer, as shown in the left panel of Fig. [21 Thus, 
a cut on (U„,-G) varying with (G-1Z) can remove most of 
the low redshift interlopers. However, the presence of dust 
in star-forming galaxies at high redshift can redden the (G- 
TZ) colour, resulting in an overlap with the colours of lower 
redshift galaxies. Thus, these high redshift dusty galaxies 
would not be classified as LBGs. 

The right panel of Fig. [!?] s hows the lumi nosity function 
of model galaxies selected using lSteidel et all magnitude and 
colour cut criteria, compared with that for all model galax- 
ies at z = 3.2. This figure shows that the lumin osity func- 
tion of galaxies selected following ISteidel et al.l criterion is 
complete above the rest-frame UV absolute magnitude that 
corresponds to the R-band magnitude limit of the observa- 
tions. Note that the observational data in the right panel of 
Fig. [2] goes deeper than the original work from lSteidel et all 

The inset in the right panel in Fig. [2] shows the pre- 
dicted redshift distribution of galaxies with 1Z < 25.5 and 
G < 27.29 and colours within the region proposed by 
ISteidel et al.l to select galaxies at 3 < z < 3.5. The model 
galaxies selected in this way have redshifts in the range 
3 < z < 4, in agreement with the selection proposed by 
ISteidel et all 



4 THE PREDICTED DUST ATTENUATION 

The presence of even small amounts of dust in high redshift 
galaxies can have a strong effect on the observ ed UV light. 
With the same model as the one used here, lLacev et al.l 
(|201ll ) predicted that dust attenuation in the far-UV (~ 
1500A) is typically 2 magnitudes for galaxies dominating 
the UV luminosity density, even out to redshift 10. 
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M ab (1500 A)-5log(h) 

Figure 3. Predicted median dust mass as a function of the atten- 
uated absolute rest-frame UV magnitude (~ 1500A) for galaxies 
at different redshifts, as indicated in the legend. For clarity, the 
error bars presenting the 10 and 90 percentiles are only shown for 
two redshifts. 

Below we present the predicted median dust masses, 
optical depths and attenuations for galaxies at 3 < z < 10. 



- z = 3.0 
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Figure 4. The predicted median difference between the attenu- 
ated and the intrinsic rest-frame UV magnitude as a function of 
the central extinction optical depth in the V-band, Ty , for galax- 
ies with Mas(1500A) - 51og/i < -17.8, plotted at the redshifts 
indicated in the legend. For clarity, the error bars representing 
the 10 and 90 percentiles are only shown for two redshifts. The 
bottom panel shows the attenuation at z = 3 for galaxies binned in 
three inclination ranges: 15° < i < 20° (solid line), 45° < i < 50° 
(dotted line, range bracketing the median value of the sample), 
70° < i < 75° (dashed line). 



4.1 The predicted dust masses 

Following iGranato et~aH |2000ri . who used a physical dust 
grain model, we assume that the dust to gas ratio in clouds 
and the diffuse ISM is proportional to the gas metallicity, 
Z coM , with a value of 1/110 for Z = Z e = 0.02. Thus, the 
total dust mass in a galaxy, Md us t, can be calculated as 
follows: 

Md - = TTo XHM - d (lof)' < 4 > 

where M co id is the total cold gas mass in the galaxy and 
Xh = 0.74 is the mean hydrogen mass fraction in the Uni- 
verse. 

Fig. [3] shows the predicted median dust mass (Eq[4]) as 
a function of absolute rest-frame UV magnitude for galaxies 
at redshifts from 3 to 10. There is a clear trend for brighter 
galaxies to have higher dust masses. This trend is a direct 
consequence of the predicted relation between the far UV 
magnitude and gas mass, combined with the fact that the 
predicted metallicity of co ld gas appears to be weakly de- 
pendant on luminosity (see lLacev et al.ll201ll ). It is also clear 
from Fig. [3] that the model predicts a decrease in dust con- 
tent with increasing redshift. At a given luminosity, we find 
a decrease in dust mass of a factor of ~5 from z = 3 to 
z = 10 for galaxies with -20 < Mab (1500A)-51og/i < -18. 

Observations of galaxies in the wavelength range from 
the far IR to the sub- mm have allowed direct estimates of the 



dust mass in some galaxies above z ~ 4. Dust masses in ex- 
cess of ~ 1O 8 /i -1 M0 have been found in galaxies with stellar 
masses estimated to be between 10 10 to 10 11 h~ 1 Mp ) at red- 
shifts 4 < z < 5 and in QSQs up to z ~ 6.5 dMcMahon et al. 
19941 : iDunlop et al.l [1994; iRobson et alj 12004 iDwek et ail 
20071 : iMichalowski et a"fll2010t ). However, other galaxies at 
simil ar redshifts do no t present any evidence of containing 
dust (|Zafar et al . 2010). For galaxies within the range of stel- 
lar masses 10 10 to 10 11 h- 1 M Q , we find median dust masses 
around ~ 10 8 h~ 1 Mq, in agreement with those observations 
which detect dust in high redshift galaxies. 



4.2 The predicted optical depth 

As described in ij2.ll the dust attenuation of starlight is di- 
rectly related to the extinction central optical depth in the 
V-band of a galaxy, tv ■ Fig. [4] presents the median pre- 
dicted attenuation in the UV as a function of tv for galax- 
ies with Mab(1500A) — 51og/i < —17.8 (corresponding to a 
te nth of the character istic luminosity at z — 3 as reported 
bv lSteidel et al.lll999h . Despite the scatter, there is a clear 
trend for galaxies with higher extinction optical depths to 
present higher UV attenuations. This trend is found to be 
independent of redshift, in the studied range 3 < z < 10. 
We find the median optical depth of galaxies brighter than 
Mab(1500A) - 51og/i < -17.8 to be between 2 and 6. These 
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Figure 5. The median attenuation in magnitudes at 1500A as 
a function of intrinsic (left panel) and attenuated (right panel) 
absolute UV (1500A) magnitude for model galaxies at z ~ 3 
(black), 6 (green) and 8.2 (red). The error bars show the 10 and 
90 percentiles for the distribution. 



median values depend strongly on the magnitude cut applied 
to select the galaxies. We also find that brighter galaxies 
present higher extinction optical depths, though this trend 
has a large scatter. 



4.3 The predicted median UV attenuation 

The attenuation of starlight depends not only on the ex- 
tinction central optical depth but also on the inclination of 
the galaxy, the distribution of the dust and the fraction of 
dust in molecular clouds. The lower panel in Fig. 0] shows 
that when galaxies are selected within small ranges of in- 
clination, the attenuation-optical depth relation shows little 
scatter, with the attenuation increasing as the inclination 
increases towards the edge-on case. 

Fig. \S\ shows the median attenuation of model galaxies 
at z = 3, 6 and 8.2, as a function of both the intrinsic and at- 
tenuated UV absolute magnitude. In both cases a trend can 
be seen for brighter galaxies to be more attenuated, though 
this is dominated by scatter, especially for attenuated UV 
magnitudes. The large scatter found suggests that, in gen- 
eral, the attenuation is not related in a simple way with the 
UV magnitude. 

Fig.[6]is similar to Fig.[5]but for the subsample of LBGs 
at z = 3, 6 and 8.2, selected by applying (to t he attenuated 
colour s) the magnitude and c olour selection of Steidel et al.l 
|l995l ). iBouwens et~aH (|2012l ) and lLorenzoni et all (|201ll ) 
respectively. Fig. [S] shows that the range of absolute magni- 
tude covered by these galaxy selection s is different. As was 
shown in Fig. [51 the LBGs selected bv lSteidel et alj criteria 
are bright galaxies, w ith Mab(1500A )-51ok/i < -20. The 
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selections from both IBouwens et al.l and lLorenzoni et alj 
make use of the deep HST photometry, sampling fainter 



M -5loq(h) M -5loq(h) 

intrinsic attenuated 

Figure 6. Figure similar to Fig. \5\ but for the subsample of 
model galaxies which pass various LBG colour and ma gnitude 
selections at z ~ 3 (selection from [pteidel et al.l 1 19951 . black) , 
6 (se lection from |Bouwe ns et al.l l 2012l . green) and 8.2 (selection 
from lLorenzoni et alj|201ll . red). For comparison, the median at- 
tenuation for all galaxies at z = 3, shown in Fig. [5] is included 
here as a grey dashed line. The error bars show the 10 and 90 
percentiles of the distribution. 



galaxies. In the next section we will explore these selections 
further, comparing them with the model predictions. 

Fig. [6] shows clearly that intrinsically brighter LBGs 
tend to be more attenuated at all redshifts. However, when 
plotting against the attenuated UV luminosity, the trend is 
weakened, being dominated by the scatter. A direct estima- 
tion of the attenuation in high redshift galaxies spanning 
the range from to 5, similar to that predi cted, has recently 
been made using results from Herschel dBurgarella et alj 
l201ll ; iBuat et al ll2012l ; iLee et al.ll2012l ; iReddv et al.ll2012l ). 

The right panel in Fig. [B] remains virtually unchanged 
if we only apply the cuts in apparent magnitude from the 
observational surveys (see the Appendix for the specific val- 
ues). This suggests that the LBG colour selection does not 
lead to an appreciable bias in average UV attenuation of 
the sample. In fact, the right panel in Fig. [6]shows the com- 
plete overlap at z = 3 between the median attenuation of all 
galaxies and that for LBGs, at least within the magnitude 
range where the latter are selected. 



5 THE PREDICTED UV COLOURS OF LBGS 

We start our study of the predicted rest-frame UV colours 
of LBGs by comparing the locus occupied by model galaxies 
with the colour region designed by observers to implement 
the Lyman drop-out technique. Usually, these colour-colour 
regions are constructed by following the colour evolution of 
model starburst galaxies with a range of negative UV con- 
tinuum slopes, predicted using a stellar population synthesis 
(SPS) model. Figs.[2l[7]and[8]show as grey areas some exam- 
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Figure 7. The left and central colu mns show the predicted density c ontours for model gala xies at the redshifts indicated in the legend, 
in the colour-colour space defined bv lBouwens et al ] (2012ft (left) and lWilkins et alj (2011bft (centre) to select LBGs at z = 5 (top row), 
2 = 6 (middle row) and z = 7 (bottom row) [grey areas in each plot] . The density contours enclose 80 per cent of galaxies brighter than 
the limits set by the different observational studies (see the Appendix). The solid/dotted contours correspond to colours with/with out 
considering the dust attenuation. The right column shows as black lines the predicted completeness and luminosity function of all galaxies 
at these redshifts. The blue /red lines correspond to the predicted comple teness and luminosity function for the subsample obtained using 
the selection criteria from lWilkins et al.1 ( 2011b) /B ouwens et alj |2012j), the predicted redshift distributions of these subsamples are 
show n in the insets of the right column. For com parison, in the right column we also show in green t he observed luminosit y function 
from: IBouwens et all i2007T ) (squ ares, 1600A) and iMcLure et al.1 i2009l ) (circles, 1500A) at z = 5 (topi; iBouwens et all j2007ll (s q uares , 
1350A1 and iMcLure et al.1 l|2009h (circles, 1500A) at z = 6 (middle): IMcLure et al] (201Cft (squares, 1500A) and lOuchi et al] |2009ft 
(circles, 1500A) at 2 = 7 (bottom). 
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Figure 8. Top row: The predicted densit y cont ours for mode l galaxies at the redshi fts indicated in the le gend, in the colour-colour 
space defined by the iBouwens et al, | (|2012| S (left) lLorenzoni et alj J201lh (centre) and IBouwens et alj J201lh (right) to select LBGs at 
z = 4, z = 8.2 and z = 10, respectively (grey areas in each plot). The density contours enclose 80 per cent of galaxies brighter than 
the limits set by the different observational studies (see the Appendix). Bottom row: The black lines show the predicted completeness 
and luminosity function of all galaxies at z = 4 (left), z = 8.2 (centre) and z = 10 (right). The red lines correspond t o the predicted 
compl eteness and luminosity functi on at the sam e redshifts for the sub sample obtained using the selection criteria from IBouwens et alj 
(2012) (left). lLorenzoni et alj 1I2OIII ) (centre) and IBouwens et alj l l201lf) (right). The predicted redshift distributions of these subsamples 
are shown in the insets (note that the y-axis units in th e left panel are No,(z)/deg 2 /1 7 ). For compa r ison, we also show in green the 
obse rved luminosity funct ion from: lYoshida et al. I J2006h (z ~ 4, squares, 1500A, left), IBouwens et alj ||2007| ) (2 ~ 4, triangles, 1600A, 
left), lMcLure et all j201Cl ) (z ~ 8.2, 1500A, centre) and Bouwens et alj (20 111) (z ~ 10, 1600A, right). 



pies of such colour selection region^. In most cases, the ob- 
servational data points are distributed covering most of the 
drop-out region. In these figures, we compare the observa- 
tionally designed drop-out regions with the density contours 
for model galaxies brighter than the limiting magnitudes in 
the corresponding observational studies (see the Appendix 
for a list of the selection criteria used in different studies). 
As we can see in Figs. [2] El and [H] the predicted UV colours 
of bright galaxies at 2.5 < 2 < 10 are consistent with the 
colour selection regions designed observationally. 

In each of the colour-colour panels in Figs. [2] [7] and [H 
the location of model galaxies at a range of redshifts, beyond 
the target one, are shown to assesses whether or not there 
are interlopers. In these figures we also present the predicted 
redshift distribution of the galaxies selected following the 



3 HST bands used here are F435W, F606W, F775W, F850LP, 
F105W, F125W, F160W (hereafter B435, V 6 06, 1775, z 8 50, Y105, 
J125, Hiao, respectively). The transmission curves of these filters 
can be seen in Fig. 1101 



magnitude and colour cuts from the observational studies 
(see the Appendix). These predicted redshift distributions 
agree with the target redshifts of each selection. 

As can be seen in Figs. O and [S] the predicted lu- 
minosity functions of those galaxies selected following the 
application of the magnitude and colour cuts described in 
the different observational studies are in very good agree- 
ment wit h observations up t o 2 ~ 10. This is remarkable, 
since the lBaugh et al.l J2005I) model was de signed consider- 
ing only LBGs at 2 = 3. iLacev et all (|201ll ) presented a de- 
tailed study of the physical processes affecting the predicted 
evolution of the UV luminosity function and a comparison 
with observational results. 

Fig. [7] shows that, with the HST filter set, it is par- 
ticularly difficult to select a complete sample of galaxies 
at z = 5 whilst at the same time attaining a minimum 
fra ction of low redshift interlopers. As previously reported 
by IBouwens et all (|2012l ). the expected locus of galaxies at 
2 = 5 in the colour-colour plot designed for the drop-out 
technique partially overlaps that for galaxies at lower red- 
shifts. Thus, in order to minimise the fraction of low red- 
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shift interlopers, the completeness of the high-z selection is 
reduced. At z = 5 we pre d ict tha t the drop-out selection 
proposed by IWilkins et al.l (|2011bh recovers ~ 70 per cent 
of all galaxies with -21.5 < Mab(1500A)< -18.5, and that 
from IBouwens et all (|2012h ~ 85 per cent. Figs. H [7] and [8] 
show that the predicted completeness of the drop-out selec- 
tion recovers more than 90 per cent of the galaxies brighter 
than the absolute magnitude implied by the apparent mag- 
nitude observational limits at z = 3, 4, 6, 7, 8 and 10. Simi- 
lar completeness are also found for the predicted galaxies at 
z = 2.5. 

Next we explore how the predicted UV colours change 
when different model parameters are varied. In particular, 
we explore the effect of changing the stellar population syn- 
thesis (SPS) model and the attenuation by both the IGM 
and dust. 



5.1 The effect of changing the SPS model 

The iBaugh et al.l (|2005h model uses by defau l t the simple 
stellar populations generated bv lBressan et a l. (1998), using 
the Padova 1994 stellar evolutionary t rackjfl and the model 
stellar atmospheres from lKuruczl (|l99# l. We have compared 
the UV colours predicted using the default SPS model wit h 
those from PE GASE .2 jFioc fc Rocca-Volmerangel Il997t l. 
iMarastonl J2005l~) andlConrov et al.l (120091 '). The SPSs models 
from both IMarastonl and Conrov et al. use stellar evolution- 
ary tracks including a model for the Thermally-Pulsating 
Asymptotic Giant Branch phase of stellar evolution that 
can reproduce observations better than pr evious attempts 
to model this stellar phase. IConrov et alj also included a 
post-AGB phase and an extended blue horizontal branch 
for stellar evolution, matching UV observations of old star 
clusters. 

The effect that stars in the Thermally-Pulsating 
Asymptotic Giant Branch phase have on the UV colours 
of a galaxy is negligible, since these stars are mainly bright 
in the near IR. Though bright in the UV, stars only reach a 
post-AGB phase at ages much older than those of the LBGs. 
Thus, little variation is expected in the UV colours predicted 
on using different SPS models. Indeed, we find that the UV 
colours of LBGs predicted using the different SPS models 
are consistent with each other. 



5.2 The effect of the IGM on the UV colours 



As described in WI?2\ the IGM attenuates the UV light, 
mainly at wavelengths A < 1216A, and has a larger effect at 
higher redshifts. Fig. [9] shows as density contours the pre- 
dicted rest-frame UV colours for galaxies at z — 4 with and 
without taking into account the effect of the IGM atten- 
uation. At z > 3, the colours sampling the Lyman-break 
change by more than 0.5 magnitudes when the IGM attenu- 
ation is included. In Fig. [9] we can see such change at z = 4. 
It is clear the need to take into account the effect of the 



4 http:/ /stev. oapd.inaf.it/cgi-bin/cmd, 
http://plciadi.pd.astro.it/ 

b These are close to the isoc hrones and stellar atmospheres used 
bv lBruzual fc Chariot! j200^ . 
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Figure 9. The (B435 -V 6 06)ab versus (V6 06-Z85o)ab colour- 
colour plot proposed bv lBouwens et al for selecting galax- 
ies at z ~ 4 (grey region) together with the predicted density 
contours for galaxies at z = 4. The density contours enclose 80 
per cent of galaxies with V606 <30.1 and zsso <29.4. The inset 
shows the predicted luminosity function for galaxies selected at 
2 = 4 with Vgoe <30.1, zgso <29.4 and the colour cuts shown by 
the grey region in the colour-colour plot. In both cases, the black 
lines correspond to the predictions without considering the IGM 
attenu ation and the b l ue/red lines to those predictions made us- 
ing the lMadaul il995r i / IMeiksinl (|2006|) prescriptions for the IGM 
attenuation. Note that the predicted luminosity functions using 
cither IGM prescription are the same at this redshift. For com- 
parison, we show in green the observed lu minosity function from 
lYoshida et"al] J2006I) (squares, 1500A) and lBouwens et al.l i2007r i 
(triangles, 1600A). 



IGM attenuation when studying the rest frame UV colours 
of high redshift galaxies. 

Fig. [9] also shows the dis tributio n of mo d el ga laxies 
when either the iMadaul (| 19951 ) or the IMeiksinl J200B) pre- 
scription is used to model the IGM attenua tion. In 32.21 we 
pointed out that the pre scription of IMadaul shows a smaller 
transmission than that of lMeiksfnl . This implies that the pre- 

i are exp e 
j IMadaul p 

tion. Fig. [5] illustrates this difference for galaxies at z = 4. 
Nevertheless, the difference in the predicted rest-frame UV 
colours is found to be small. In fact, we do not find a signifi- 
cant difference in the luminosity function of LBGs modelled 
using either IGM attenuation prescription for galaxies at 
2.5 < z < 10 (see the inset in Fig. [9] for the case of z = 4 
and LBGs selected following IBouwens et al.ll2012h . 

In our model, reionisation is assumed to be instanta- 
neous. Indeed, analyses of the line of sight to quasars at 
2 > 6 are consistent with a uniformly ionised IGM at z ~ 6 



dieted colours sampling the Lyman-break are exp ected to 
be slightly redder for the model using the IMadaul prescrip- 
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(e.g. ISchroeder et all 120121 ). At higher redshifts, numerical 
simulations of reionisation suggest a gradual process, im- 
plying t he existence of ionised bubbles in a mainly neutral 
medium |Benson et al.ll200ll ; liTiev et al]|2006l ). The fact that 
the predicted luminosity function including the effect of the 
IGM matches the observations suggests that a minimal frac- 
tion of galaxies would not be detected due to the possible 
patchiness of the IGM at very high redshifts. 

5.3 Dust treatment and the UV colours 

In Figs. [7] and [8] we compare the predicted distributions of 
UV colours with and without considering the effect of dust. 
As can be seen there, in general, the presence of dust in- 
creases the spread of UV colours. At z = 5, the inclusion 
of dust clearly makes galaxies redd er in the colou r sam - 
plin g the Lyman-break fo r both the iBouwens et al.l (|2012l ) 
and lwilkins et all (|2011bh colour selections. At this redshift 
most galaxies would be missed by these two drop-out tech- 
niques if intrinsic colours were going to be used. As men- 
tioned before, using the HST bands, the selection of LBGs 
at z = 5 needs to be done using stringent colour cuts, in or- 
der to minimise low redshift interlopers. The effect of dust 
on the UV colours at other redshifts, z ^ 5, is not so strong 
and approximately the same completeness is recovered using 
the predicted galaxy colours with and without considering 
the effect of dust. 

Next we explore the effects that changing some of the 
model parameters controlling the dust treatment have over 
the predicted rest-frame UV colours. 
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Figure 10. Top panel: The optical depth as a function of wave- 
length for 8 galaxies at z = 4 with different UV attenuations, ob- 
tained using 42 narrow top-hat filters starting with cither the MW 
extinction curve (solid lines) or the SMC curve (dashed lines). 
Bottom panel: The optical depth normalised at 1.5/im for those 
galaxies in the top panel (each galaxy is presented as a line of 
the same colour in both panels). For reference, also shown are 
the response functions of the HST filters used to select LBGs at 
different redshifts (lower set of grey curves). 



5.3.1 Changing the extinction law 

The extinction curve of a galaxy is determined by the 
composition and size distribution of the galactic dust, 
which in turn depends directly on the star formation 
history of a galaxy. At high redshift it is particu- 
larly unclear w hich mechanism produced t he bulk of the 
obser v ed dust dBianchi fc Schneiderl 20071; Valiante et al.l 



20091; iMichalowski et all |2010| ; iDwek fc Cherchnefj l201ll ; 
Valiante et al.ll201ll ) and thus which extinction curve best 
describes high redshift galaxies. Here, the dust attenuation 
calculatio n is based on the re sults of the radiative transfer 
model of iFerrara et al.l (1 19991 ), which assumes a particular 
extinction curve as an input. 

Fig. [10] shows for 8 galaxies, the effective optical depth 
as a function of wavelength, obtained using 42 narrow, 50A, 
top-hat filters starting with either the MW extinction curve 
(the default one in the model) or the SMC extinction curve. 
These 8 galaxies are at z = 4, they are bright, zsso > 29.4, 
they have inclinations between 55 and 65 degrees and they 
have been selected as they exhibit different attenuations in 
the ^850 band. We can see from this figure that the main fea- 
tures of both the MW and SMC extinction curves, namely 
the bump at 2175A for the MW curve and the stronger UV 
extinction at smaller wavelengths in the SMC curve, are still 
visible after the radiative transfer calculation, which takes 
into account the distribution of dust in the galaxies and the 
viewing angle. As we can see in Fig. 1101 the characteristics of 
the assumed extinction curve affect the inferred dust atten- 
uation. This result stresses the relevance of understanding 



better the characteristics of dust distribution when inter- 
preting observational results. 

The lower panel in Fig. [10] shows the effective optical 
depth normalised to that at 1.5/im, when starting with a 
MW extinction curve. This panel shows clearly how the fea- 
tures in the attenuation curves get smoother for higher at- 
tenuations. The same is seen if we start with a SMC extinc- 
tion curve for dust in both the diffuse phase and in molecular 
clouds. 

Fig. [11] shows the distribu tion of model galaxies in the 
colour-colour plot proposed by IBouwens et al.l (|2012l ) to se- 
lect galaxies at z ~ 6 (grey area) when different input ex- 
tinction curves are used. The model galaxies populate the 
region designed to select galaxies at z ~ 6, independently of 
the input extinction curve. However, the median UV contin- 
uum slope is redder if we start with a SMC extinction law. 
We will return to this point in section i}6] 

The inset in Fig. [TT] shows th at the luminosity fun ction 
of LBGs at z=6 (selected following IBouwens et al.l(2012l ) pre- 
dicted from the model using the SMC extinction law is about 
half a magnitude fainter than that using the MW extinction 
law. The luminosity function of all galaxies at z — 6 starting 
the calculations with either the SMC or the MW extinction 
curves still show this difference and, thus, it cannot be due 
to the colour cuts applied in the LBG selection. This differ- 
ence is due to the fact that the SMC curve, as shown in Fig. 
[JJ produces a larger attenuation in the rest-frame UV range 
than is obtained fro m the MW curve. 

As discussed in lLacev et all (|201ll ), the UV luminosity 
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Figure 11. The (i775- Z85o)AB versus (z85 -Ji 25 )a.b colour- 
colour plot proposed by iBouwens et al.l || 2012T ) to select galax- 
ies at z ~ 6 (grey region) with the density contours for model 
galaxies at z = 6. The density contours enclose 80 per cent of 
galaxies with zgso < 29.4 and J125 < 29.9. The inset shows the 
predicted luminosity function for galaxies selected at z = 6 with 
Z850 < 29.4 and J125 < 29.9 and the colour cuts shown as a grey 
region in the colour-colour plot. In both cases, the black lines 
corresponds to the default case of using the MW extinction law, 
the red to using the SMC extinction law and the blue lines cor- 
responds to predictions made using the SMC extinction law and 
adopting / c ioud = 0-7 (where the default value is / c ioud = 0.25). 
For comparison, we show in green the obs ervational data from 
IBouwens et al l (|2007f l (squares, 1350A) and lMcLure et alJ l|2009T ) 
(circles, 1500A). 

function is very sensitive to the fraction of dust in molecular 
clouds, / c ioud, and the time that is assumed for stars to mi- 
grate out of clouds, icsc- Increasing /cioud or decreasing i eS c 
has the effect of reducing the dust attenuation. By default, it 
is assumed that 25 per cent of the dust is locked up in molec- 
ular clouds and that new stars will migrate out of these after 
lMyr. From MW observa tions, this escape time is fo und to 
be between 1 and 3 Myr (jLeisawitz fc Hauserlll988l ). Thus, 
further reducing t csc might be unrealistic. However, increas- 
ing the amount of dust in clouds from 25 per cent could be 
reasonable for some galaxies, in p articular, for star form- 
ing galaxies at high redshift (e.g. iDaddi et a.1.1 l201oj ) . The 
inset of Fig. 1111 shows that assuming / c ioud = 0.7 increases 
the amplitude of the UV luminosity function of LBGs mod- 
elled with the SMC extinction law, so it gets closer to the 
prediction obtained using the MW extinction law and there- 
fore to the observational data. However, the predicted galax- 
ies still populate a similar region in the colour-colour plane 
shown in Fig. [11] and we find a UV continuum slope con- 
sistent with that obtained starting with the default value of 
/cloud = 0.25. 

We have tested that changing either / c i ou d or t csc has 
a negligible impact on the predicted UV colours of model 



galaxies. However, the luminosity function depends on these 
parameters and thus, the number of selected LBGs will be 
different for different sets of parameters. This implies that 
UV colours cannot constrain the dust par ameters, /cioud 
and t e sc . A similar conclusion was reached bv lLo Faro et all 
(|2009h in their study of model LBGs at 4 < z < 6. 



5.3.2 Changing the spatial distribution of dust 

The way dust is distributed in a galaxy with respect to stars 
affects the resulting attenuation and thus, it might affect the 
predicted UV colours for a given galaxy. By default in GAL- 
FORM it is assumed that the scale height of the dust distri- 
bution follows exactly that of the stars, h Zt dust/h z , s tars = T 
In order to explore how this affects the predicted UV colours 
and luminosity functions of LBGs at 3 < z < 8.2, we have 
made predictions for two other cases: when dust is more con- 
centrated than stars, h Zt dust/h Z: stars — 0.4, and when it is 
more extended, h Z}duat /h z>s tars = 2.5. 

The UV luminosity function is affected by the change 
in the relative scale heights of dust and stars. When dust is 
more concentrated, it affects the light of a smaller fraction 
of stars and thus the galaxy presents a smaller attenuation 
compared to the default case of having h Zt d U st/h Zta tars = 1. 
When the dust has a larger height distribution than the 
stars, all of the starlight will be affected by dust, as happens 
in the case of where h Zt dust/h ZtS t ar s = 1, and thus, the two 
luminosity functions are very close. 

Despite the change in the UV luminosity function, the 
model galaxies still populate approximately the same region 
in the UV colour-colour planes. 



6 THE UV CONTINUUM SLOPE 

The rest-frame UV continuum follows approximately a 
power law, fx oc X 13 , and thus , it can be characterised by 
the slope /3 i|Meurer et"ai]|l999l fl. The UV continuum slope 
is often considered to be directly relate d to the amount of 
dust in a galaxy. In lWilkins et all (|2012h we studied the pre- 
dicted intrinsic (without attenuation) UV continuum slope, 
measured from a single colour. We found that the distribu- 
tion of intrinsic UV continuum slopes is affected by the star 
formation and metal enrichment histories of galaxies and by 
the choice of the IMF. These dependencies limit the accu- 
racy of the estimation of dust content in galaxies from the 
observations of UV continuum slopes. We also found that 
the intrinsic UV continuum slope is correlated with differ- 
ent galaxy properties, including redshift, suggesting that the 
attenuation by dust inferred from observations may be mis- 
estimated, if these trends are not taken into account. 

In most observational studies until now the UV contin- 
uum slope, fj, has been measured from a single UV colour. 
Although this method is hindered by the photometric errors 
on the two bands used, a fit to the spectral energy distri- 
bution of g alaxies at high redshifts is only now be coming 
possible (|Finkelstein et al.ll201ll ; iDunlop et al.||2012T ). Thus, 



6 Note that a /3 = — 2 corresponds to a galaxy with a flat spec- 
trum in terms of /„ and, thus a zero UV colour in AB magnitudes. 
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Figure 12. The UV continuum slope of model galaxies at z = 5 
(solid lines) and 2 = 7 (dotted lines) measured with two top-hat 
filters, NUV (~ 1500A) and FUV (~ 2200A), as a function of the 
rest-frame UV absolute magnitude. The black line corresponds 
to the intrinsic UV continuum slope and the blue/red to that 
obtained using the MW/SMC extinction curve as the input for 
the dust attenuation calculation. 



in order to facilitate a comparison with different observa- 
tional studies, we calculate ft using a single UV colour. 

We have studied the impact on the UV continuum slope 
of modifying all the model parameters discussed in ij5] We 
find that the main sensitivity is to the input extinction 
curve, which we discuss further below. As expected from the 
previous sections, changing the SPS model does not affect 
the predicted UV continuum slope. We have previously seen 
that the IGM attenuation modifies the UV colours, since this 
effect acts to enhance the strenght of Lyman-break. How- 
ever, the effect that the IGM attenuation has on the UV 
continuum slope is negligible. We have studied how the UV 
continuum slope is modified when the scale height of dust is 
different from that of stars. As we saw in H5,3.2l such a change 
mainly modifies the UV luminosity function. Accordingly, if 
no magnitude cut is applied, we find no variation in the 
UV continuum slope when we go from h Zt dust/h z , s tars = 1 
to h Zi dust/h z , stars = 0.4 or to h z , dust I h z , stars = 2.5, in the 
studied redshift range. Similar results are found when chang- 
ing the parameters controlling the fraction of dust in clouds 
and the migration time out of clouds for young stars. 



6.1 The change of the UV continuum slope with 
luminosity 

Fig. [12] shows the median UV continuum slope for galaxies 
at z = 5 and z = 7 '. The UV continuum slope, ft, has been 
obtained using: 



ft = 2.47(FUV - NUV) 



(5) 



where (FUV-NUV) is the rest-frame UV colour calculated 
with two top hat filters covering the wavelength range: 
1300 < \ F uv{k)< 1700 and 1800 < W (A)< 260(Q. 

As reported in lWilkins et al.l (|2012T ), our galaxy forma- 
tion model predicts an intrinsic UV continuum slope that 
is practically independent of the UV luminosity of galaxies. 
We find that the predicted UV continuum slope varies by 
less than 0.2 magnitudes for —22 < A/tjv — 51og/i < —16. 
We find an intrinsic UV continuum slope of about —2.5 for 
galaxies 5 < z < 7. Similar values have been found in hydro- 
dynamical simulations developed to study LBGs at z > 5 
i|Forero-Romero et alJliuTol : iDaval fc Ferrarall2012r i. 

Fig. [T2] also shows the predicted UV continuum slope 
once the effect of dust is taken into account, starting from 
either the MW or the SMC extinction curves. We find that 
the UV continuum slope is very sensitive to the extinction 
curve used as input for the dust attenuation calculations. For 
the chosen set of filters, we find that, as expected, galaxies 
get redder when the dust attenuation is calculated starting 
with the SMC extinction curve. However, as shown in Fig. 
1121 model galaxies get bluer when the dust attenuation is 
calculated starting with the MW extinction curve. 

This behaviour is directly related to the 2175A bump 
in the MW extinction curve shown in Fig. [1] that is still 
present once the geometry of the dust is taken into ac- 
count, as seen for the 8 example galaxies in Fig. 1101 Thus, 
this behaviour will depend on the particular set of filters 
chosen to measure the galaxy colour and it is expected 
to be reduced if a SED fitting is used to ca l culate the 
UV continuum slope. B oth iForero- Romero et al.l (|2010l ) and 
IDaval &c FerrarJ (|2012T ) estimated the UV continuum slope 
by SED fitting their model galaxies, finding very good agree- 
ment with the observed values. These two studies estimated 
the dust att enuation assuming the simple approximation of 
a dust slab. IDaval fc Ferraral (|2012h assumed that the dust 
present in galaxies at z > 6 is produced solely by Type-II 
SNe. The extinction curve of dust produced solely by SNe 
pres ents a bump at ~ 2500A a nd a very steep slope in the 
UV |Bianchi fc Schneiderll2007l ). These characteristics could 
also imply that galaxies with dust solely produced by SNe 
will appear to having bluer UV continuum slopes if they 
were estimated using a single colour. 



6.1.1 Comparison with observations 

In order to directly test our predictions against observations, 
we calculate ft using the exact filters and equations proposed 
in different studies that use a single UV colour to estimate 
the slope of the UV continuum. Fig. H3l shows the predicted 
intrinsic UV continuum slopes calculated following different 
observational studies for galaxies at z ~ 5, 6 and 7, using the 
magnitude and colour cuts summarised in the Appendix. 

At z = 6, Fig. [13] illustrates how, when starting from 
the MW extinction curve, model galaxies can look redder, 
the same or even bluer than their intrinsic colours, depend- 
ing on the choice of filters. This is directly related to the 
2 175 A bump being sampled or not by the two chosen filters. 
This extreme sensitivity of ft to the choice of colour used in 



7 Note that these two filters cover the same wavelength range as 
the GALEX niters. 
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Figure 13. The predicted median intrinsic UV continuum slope (black lines) and the attenuated slope, starting either with the MW 
(blue lines) or the SMC extinction curve (red lines), as a function of absolute attenuated UV magnitude, for galaxies at z ~ 5 (top), 
z ~ 6 (middle row), z ~ 7 (bottom), selected following the magnitude and colour cuts proposed in different observational studies 
(summarised in the Appendix). The different columns correspond to the UV continuum slope calculated following the method adopted in 
the observational studies indicated on the top of each plot. The error bars show the 10 and 90 percentile range of the model pre dictions. 
The green ci rcles present, for com parison, the observational data from lBouwens et al.l ||2012T ) (left column) , IWilkins et alj J201 lbl ') (middle 
column) and lMcLure et all l|201lh (right column, who selected LBGs based on photometric redshift and not UV colour cuts). In the latter 
case we also show the observations for individual galaxies in light green. The mean observational points are shown with their associated 



the es timation reinforces the results from iFinkelstein et all 
l|201l| ). who showed how weakly the UV continuum slope is 
constrained with a single colour compared with SED fitting 
using at least 5 bands. 

The right column of Fig. \T3\ shows the UV continuum 
slopes obtained for individual galaxies in the iMcLure et al.l 
(|2011T ) sample, together with the mean val ues per magni- 
tude b in, as shown for the s amples from both lBouwens et alJ 
(|2012h and IWilkins et all (|2011bl ) in the other columns of 
Fig. [13] The scatter when considering individual galaxies is 



large enough to cover the UV continuum slopes predicted 
starti ng with both the MW and SMC extinction laws. How- 
ever, IMcLure et ail (|20Tll ) presented a detailed analysis of 
the problems that affect individual measurements and con- 
cluded that averaged UV continuum slopes are more reliable. 

Fig. [13] shows that, for some filter combinations, the 
predicted /3 including the effect of dust agree with the ob- 
servations. However, in most cases, the observed mean val- 
ues are closer to the predicted intrinsic values than to those 
including the effect of dust. 
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By doing the exercise of smoothing out th e 2 175 A bump 
in the attenuation curves calculated with the iFerrara et all 
1 19991 ) model (which are still a function of inclination, opti- 
cal depth, etc.), we obtain UV continuum slopes about 0.2 
magnitudes redder than the intrinsic ones and which are 
therefore closer to, and in some cases consistent with, the 
observed values. We have also used the attenuation at 1500 A 
obtained in the above calculation to normalise the applica- 
tion of the Calzetti law to the model galaxies to derive an- 
other illustrative estimate of the attenuation. We find that 
the attenuation curves from this procedure are greyer than 
the Calzetti law i.e. there is less wavelength dependence. 
By using the Calzetti law normalised in this way, we ob- 
tain attenuated UV continuum slopes which are consistent 
with those inferred observationally in at least one magni- 
tude bin at z=5, 6 and 7. Therefore, assuming the Calzetti 
law, observations are consistent with the predicted 1500A 
attenuations. 

The predicted values of the U V continuum slope are ex- 
tremely sensitive to the calculation of the dust attenuation. 
The observational data are well within the range of predic- 
tions when considering either the MW or the SMC extinc- 
tion curves as input, again emphasising the uncertainty in 
the treatment of attenuation, and the implications that has 
when interpreting observational data. It is unclear which ex- 
tinction curve should be used at different redshifts or even at 
different luminosities. The fact that the mean observational 
data have values between the intrinsic model predictions for 
P and those obtained using the SMC extinction curve, might 
suggest that the MW extinction curve is not suitable for de- 
scribing galaxies at z > 5. Nevertheless, the use of the SMC 
extinction curve appears to be also disfavoured since model 
galaxies are too red compared with observed ones. Our re- 
sults suggest that the extinction curve that might describe 
high redshift galaxies will h ave characteris t ics be tween those 
of the MW and SMC ones. IConrov et all (|2010h found that 
the colours of galaxies at z < 1 are best described if it is as- 
sumed an attenuation curve with a UV bump strength equal 
to 80 per cent of that of the MW and a higher effective op- 
tical depth in the UV range is assumed, compared with the 
MW curve. 

One clear point of tension in Fig. [13] is the ob- 
served trend for f3 to get bluer for fainter galaxies, some- 
thing that is not seen i n the model galaxies. Nevertheless, 
iFinkelstein et all (|201ll ) have argued that this trend is in 
fact due to assuming an average redshift for the whole sam- 
ple and applying a single correction to estimate j3 as a 
function of the re st-frame UV luminosity using fixed filters. 
IFinkelstein et all found that once a more consistent defi- 
nition of the rest-frame UV luminosity is used, the trend 
of the UV continuum slope wit h UV luminosity d isap- 
pears. The observations from both lDunlop et all (|2012T ) and 
IFinkelstein et al.l (|201lh resulted in mean UV continuum 
slopes independent of the galaxy luminosity, in agreement 
with our predictions. These two observational studies se- 
lected LBGs by measuring their photometric redshifts and 
used the galaxy SED to obtain their absolute rest-frame UV 
magnitude. 

For the magnitudes measured in a particular band we do 
find that brighter galaxies are, on average, more attenuated. 
However, the model predicts rather grey attenuation and so, 
the predicted attenuations in the two bands used to compute 



q 

CN 

I 

3 



> 
Z) 
l_i_ 

sf 

CN 

II 








8 



Figure 14. The UV continuum slope of galaxies with rest- 
frame Myv — 51og/i < —20.3 (solid lines, £Jj V at z = 3) and 
M\jv — 51ogh < —17.8 (dotted lines, O.lLJjy at z = 3) measured 
with two rest-frame top-hat filters, NUV (~ 1500A) and FUV 
(~ 2200A), as a function of redshift. The black line corresponds 
to the intrinsic UV continuum slope, the blue/red to that ob- 
tained using as an input for the model the MW/SMC extinction 



P are very similar for different HST bands, giving attenuated 
colours that practically do not depend on the luminosity of 
a galaxy. 

Fig. [13] shows that extracting information about the 
dust content of a galaxy based on its UV continuum slope 
could lead to erroneous conclusions if knowledge of the size 
distribution of dust grains and the galactic inclination is not 
available. 



6.2 The change of the UV-continuum slope with 
redshift. 

Fig. [14] shows the evolution with redshift of the intrinsic 
UV continuum slope, calculated with Eq[5] together with 
that estimated starting from either a MW or a SMC extinc- 
tion curve for galaxies with Mab(1500A) — 51og/i < —20.3 or 
Mab(1500A) — 51og/i < —17.8. These two rest-frame mag- 
nitudes correspond to V uv and 0.1L£,^ galaxies at z = 3 
(|Steidel et al.ll 19991 ). 

The intrinsic UV continuum slope becomes bluer (i.e. 
P becomes more negative) with increasing redshift. For the 
intrinsic UV continuum slope, we find a change of 0.5 mag- 
nitudes from z — to z = 8, independently of the ab- 
solute magnitude cut. This change is due to an increas- 
ing population of young and massive stars at higher red- 
shifts. This result agrees with most observational studies 
ilBouwens et alJl20od ; IFinkelstein et al.ll201ll ; IWilkins et~al] 
120121; iBouwens et alj|2012h . but disagrees with the results 



from iDunlop et al.l (|2012l ) for LBGs at 5 < z < 7. The fact 
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that the intrinsic UV continuum slope for model galaxies 
varies with redshift suggests that the trend found observa- 
tionally may not entirely be driven by the evolution in the 
dust content of galaxies. 

Fig. [14] shows that the effect of using either the SMC 
or the MW extinction curve impacts the predicted UV con- 
tinuum slope in a similar way to that reported previously 
for a fixed redshift. When starting with the MW extinction 
curve we still find that the UV continuum slope becomes 
bluer with increasing redshift. However, when starting with 
the SMC extinction curve the variation of UV continuum 
slope with redshift is negligible. 



7 CONCLUSIONS 

In this paper we have continued the study started by 
lLacev et al.l |201ll ) of Lyman-break galaxies (LBG) using 
GALFORM, the semi-anal ytical mod el of g alaxy formation 
originally developed by ICole et"aTl (|2000l ). We have fol- 
lowed the galaxy formation process in the framework of the 
ACDM cosmology. We have presented predictions for the 
iBaugh et al. I (|2005h variant of GA LFORM, with the sa me in- 
put parameters as those used in lLacev et all (|201ll l This 
model was originally constructed to match the far-UV lu- 
minosity function of LBGs only at z = 3 but was shown to 
fit the observed luminosit y functions up to z — 10, the full 
range currently available (|Lacev et al.ll201lh . 

Here we have studied the predicted re st-frame UV 
colours of LBGs in the range 2.5 < 2 < 10. The lBaugh et all 
model of galaxy formation includes physical treatments of 
the hierarchical assembly of dark matter haloes, shock- 
heating and cooling of gas, star formation, feedback from su- 
pernova explosions, the photoionization of the intergalactic 
medium (IGM), galaxy mergers and chemical enrichment. 
The luminosities of galaxies are calculated from a stellar 
population synthesis (SPS) model and dust attenuation is 
then included using a self-consistent theoretical model based 
on the results of radiative transfer calculations. The far-UV 
dust attenuation is a critical component in any model for 
LBGs, since the effect of dust on UV light can be dramatic. 

We find that model galaxies have dust masses that 
decrease for fainter galaxies. Model galaxies with —22 < 
Mab(1500A) — 51og/i < —18 have median dust masses rang- 
ing from lO 8 5 fe _1 M to lO^h^M® at z = 3. We also find 
that median dust masses decrease with increasing redshift. 

We find that brighter LBGs are generally predicted to 
be more attenuated at all redshifts than fainter ones, though 
this trend is dominated by the scatter. Our results suggest 
that the LBG colour selection does not lead to an apprecia- 
ble bias in average UV attenuation of the sample. 

Our model produces galaxies with UV colours that are 
generally consistent with the colour-colour region designed 
observationally for selecting LBGs. From the model results, 
we expect the drop-out technique to recover most galaxies 
brighter than the absolute UV magnitude corresponding to 
the apparent magnitude limit of a particular observational 
survey. We find the completeness of the drop-out technique 
to be above 90 per cent at all the studied redshifts, ex- 
cept z = 5, at which we predict a completeness above 70 
per cent. These values should be considered as upper limits, 
since we have not included the effect of photometric noise, 



which most likely will push some galaxies out of the selection 
region at the target redshift. 

We have investigated the impact that different parame- 
ters from the model (SPS, IGM and dust treatment) have on 
the UV colours, finding that they are most sensitive to dust 
and, in particular, to the extinction curve assumed as an in- 
put for the radiative transfer model with which attenuation 
curves are obtained for individual galaxies. This is related 
to the fact that the characteristics of the input extinction 
curves are recognisable in many of the attenuation curves 
obtained after the p rocessing performed with the radiative 
transfer model from iFerrara et all (| 1999h . Nevertheless, we 
find that the drop-out technique that selects LBGs based on 
two UV colours is very robust. 

The predicted UV luminosity function is strongly af- 
fected by changes in the model parameters t osc , the time 
that it takes a star to escape from the molecular cloud in 
which it is born, and /cloud, the fraction of gas and dust in 
molecular. However, we find that the predicted UV colours 
are not sensitive to these two parameters. Thus, although 
the observed UV colours of LBGs are reproduced by the 
model, the UV colours do not help further to constrain the 
dust parameters t csc and / c ioud- 

We have also obtained the predicted UV continuum 
slope calculated from a single rest-frame UV colour, as is 
commonly done in observational studies. For galaxies at red- 
shifts between 5 and 7 we find an intrinsic (without any at- 
tenuation) UV continuum slope of about —2.5, even though 
some of the brightest galaxies have up to 2 magnitudes of 
extinction in the UV. 

Using the Milky Way (MW) dust extinction curve as 
an input to the radiative transfer model that calculates dust 
attenuation, the predicted UV continuum slopes are, in gen- 
eral, bluer than observations. However, we find the opposite 
trend when the Small Magellanic Cloud (SMC) extinction 
curve is used as an input. This demonstrates the strong de- 
pendence of UV colours on dust properties and the inade- 
quacy of using the UV continuum slope as a tracer of dust 
without further knowledge of the galaxy inclination or dust 
characteristics. The mean observed UV continuum slopes 
are close to the predicted intrinsic ones. Including the at- 
tenuation by dust produces UV continuum slopes that are 
either too blue, if assuming a MW extinction curve, or too 
red, if assuming an SMC extinction curve, compared with 
the observations. This suggests that neither of these extinc- 
tion curves provides an adequate description of high red- 
shift galaxies. The discrepancy between the model predic- 
tions and the observations for the slope of the attenuated 
UV continuum slope is smaller than the difference between 
the predictions obtained using either the MW or SMC ex- 
tinction curves. 

The predicted median UV continuum slope depends 
very weakly on the rest-frame UV luminosity of galaxies. 
This result agrees with the observational estimates from 
both lDunlop et all (|2012i ') and lFinkelstein et all J201ll). but 
disag r ees with other obser v ational studies llBouwens et all 
2009*; IWilkins et al. I l2011bl ; iBouwens et~all I2012T I .The pre- 
dicted UV continuum slope gets bluer with increasing 
redshift, in agreement with many observational studies 
ilBouwens et alll2009 |; iFinkelstein et al.ll201ll ; IWilkins et ail 
120121 ; IBouwens et alj|2012l ). though this trend is not clear 
when starting with the SMC extinction curve. 
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In conclusion, we find that the iBaugh et all (|2005l ) 
model predicts UV colours consistent with the colour-colour 
regions designed observationally for selecting LBGs. The 
predicted intrinsic UV continuum slopes are very close to 
those observed. However, once the (uncertain) dust attenu- 
ation is modelled, the predicted and observed UV continuum 
slope are inconsistent. The discrepancy is extremely sensi- 
tive to the choice of input dust extinction curve for the ra- 
diative transfer model that calculates the dust attenuation. 
Actually, it is unclear which dust extinction is appropriate 
for high redshift galaxies. A better knowledge of dust prop- 
erties in galaxies at high redshifts will be required to further 
constrain the models using LBGs. 
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APPENDIX 

Here we define explicitly the magnitude and colour cuts used 
to select the model LBGs at different redshifts. The mag- 
nitude cuts have been chosen to be either the magnitude 
cutoff used in the observational studies to define their LBGs 
sample or the 3a (or 5a) detection limits quoted in those 
studies for the bands sampling the UV continuum slope in 
the deepest field of the study. 

• U-dropouts, 3 < z < 3.5, ISteidel et~ail (|l995T ) 
selection: 1Z < 25.5 and G< 27.3 and 

(U„-G)> 1.5 and (U„-G)>(G-ft)+1.5 and 
(G-TZ)< 1.2 and (G-TZ)> 

• B-dropouts, 2 ~ 4, iBouwens et all (j 20121 ) selec- 
tion: Veoe < 30.1 and z 8 5o < 29.4 and 

(V 6 o6-Z85o)< 1-6 and (B 4 35-V 6 06)> 1-1 and 

(B435 — V606 ) > ( V606 — Z850 ) + 1 ■ 1 

• V-dropouts, 2 ~ 5, IBouwens et all (|2012h selec- 
tion: 1775 < 29.9 and z$ 5 o < 29.4 and 

[(Vaoe -1775) > (1775 -Z6so)+1.5 or (Vco6-i77s)> 2] and 
(i775-Z65o)< 0.8 

• V-dropouts, z ~ 5, IWilkins et al.l (|2011bT ) selec- 
tion: Z850 < 27.9 and 

(z 8 50-Y 105 )< 0.62 and (V 6 06-i77 5 )> 0.81(z 850 -Yio5)+2.14 
and (i 7 75-zg 5 o)< 0.76(z S50 -Yio5)+0.83 

• i-dropouts, 2 ~ 6, IBouwens et al.l ([201$ ) selec- 
tion: zgso < 29.4 and J125 < 29.9 and 

(i 77 5-zs5o)> 1.3 and (z 85 o-Ji25)< 0.9 

• i-dropouts, 2 ~ 6, IWilkins et~aD (|2011bh selec- 
tion: Y105 < 28.1 and (Y105-J125K 0.78 and 
(i 77 5-z 8 5o)> 0.56(Yio5-Ji2 5 )+1.27 

and (z 850 -Yio5)< 0.66(Yi 05 -Ji2 5 )+0.73 

• z-dropouts, 2 ~ 7, IBouwens et all (|20120 selec- 
tion: Y105 < 29.6 and J125 < 29.9 and 

(z 8 50-Yio 5 )> 0.7 and (Yi 05 -Ji2 5 )< 0.8 and 

(z 8 50-Yi 05 )> 1.4(Yl05-Jl2 6 )+0.42 

• z-dropouts, 2 ~ 7, IWilkins et al.l (|2011bh selec- 
tion: J125 < 28.5 and (J125— Hieo)< 0.95 and 
(Yi 05 -Ji2 5 )< 0.46(Ji25-Hi 60 )+0.89 

and (z 850 -Yio5)> 0.42(Ji2 B -Hi 6 o)+1.48 

• Y-dropouts, 2^8, lLorenzoni et all (|201ll ) 
selection: J125 < 28.5 and 

(Yio5-Ji25)> 0.9 and (J125— Hi 60 )< 1.5 and 
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(Y 105 -Jl25)> 0.73(Jl25-Hl 60 )+0.9 

• J-dropouts, z ~ 10, iBouwens et alj (|201lft selec- 
tion: Hieo < 29.8 and (Ji2 5 -Hi 6 o)> 1.2 
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